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On the kHz QPO frequency correlations in bright neutron 
star X-ray binaries 
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ABSTRACT 

We re-examine the correlation between the frequencies of upper and lower kHz quasi- 
periodic oscillations (QPO) in bright neutron-star low-mass X-ray binaries. By in- 
cluding the kHz QPO frequencies of the X-ray binary Cir X-l and two accreting 
millisecond pulsars in our sample, we show that the full sample does not support the 
class of theoretical models based on a single resonance, while models based on rela- 
tivistic precession or Alfven waves describe the data better. Moreover, we show that 
the fact that all sources follow roughly the same correlation over a finite frequency 
range creates a correlation between the linear parameters of the fits to any sub-sample. 
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1 INTRODUCTION 

Thanks to the Rossi X-Ray Timing Explorer (RXTE), 
the past decade has seen the discovery of high-frequency 
quasi-periodi c oscillations (QPO s) in neutron-star X-ray 
binaries (see Ivan der Klisl I2006L for a review) . These os- 
cillations often appear in pairs, with frequencies ranging 
from a few hundred Hz to more than 1 kHz, hence the 
name kHz QPOs. Such fast variability provides a probe 
into the accretion flow very close to the compact objects 
and is a powerful tool to observe effects of general rel- 
ativity. The frequencies of the kHz QPOs are strongly 
correlated with oth e r tim i ng and spectral features (see 
Ford fe van der Klisl 1 19981: iPsaltis. Belloni fc van der Klisl 



19991: iBelloni. Psaltis. fc van der KlisJ 

van Straaten. van der Klis fc Wiinandsl 2004 ) . 



2002; 



Several models have been proposed for these oscil- 
lations; these models are based on the identification of 
the QPO frequencies with various characteristic fr e quen- 
cies in the inner accretion flow ( Stella fc Vietri 19991; 
Osherovich fc Titarchuk|[l999l ; IZhanell2004 lLamb fc Miller! 



20031 ). However, there is still no consensus as to the origin 



of the kHz QPOs. 

QPOs in the range 30 — 450 Hz have been observed 
in a limited number of black-hole candidates; in some of 
these systems the QPOs seem to move in frequency (see 
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e.g. iHoman et al.l l200ll ). However, in the four cases when 
two QPOs are detected simultaneously in a given black-hole 
candidate, these QPOs appear at frequencies that are consis- 
tent with being always the same, with frequency ratios that 
are consistent with the ratio of two int eg er numbers, such a s 
3:2 and 5:3 (see e.g. lStrohmaverll2001al lbt lMiller et all200ll ). 
A model for these QPOs based on the resonance of charac- 
teristic frequenci es at a specific radius in the acc retion disc 
was presented by I Abramowicz fc Kluzniakl (|200ll ). 

Hb ramowicz et al. ( 2003) reported that the ratios of the 
frequencies in the kHz QPOs from the brightest low-mass 
X-ray binary (LMXB) in the sky, Sco X-l, tend to cluster 
around a value of 1.5, which they interpret as evidence for 
a 3:2 resonance also being responsible for the QPO pairs in 
neutron-star systems. 

IBelloni. Mendez fc Hornanl (|2005l . hereafter BMH05), 
have shown that the frequencies of the lower, ui ow , and 
upper, Uhigh, kHz QPOs are strongly correlated along a 
line which is inconsistent with Vhigh = l-5vi ow , but can be 
roughly approximated as a straight line not passing through 
the origin. They also showed that when this correlation is 
taken into account, the pe ak in the distribution of ratios in 
the Sco X-l data used by I Abramowicz et al.l (|2003l ) is only 
marginally significant, since the distribution of the ratio of 
two correlated quantities is completely determined by the 
distribution of one of them. BMH05 also found that the fre- 
quency distributions of lower kHz QPOs for five systems 
have statistically significant peaks (some of which did not 
occur near a 3:2 frequency ratio). However, those peaks are 
consistent with being caused by a (nearly) random walk of 
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the QPO frequency with time, together with the sampling 
of the observations. 

A math ematical approa ch to the resonance model was 
presented by Rebuscd (|200J) in order to explain the Sco X-l 
results; iRebuscol l|20fM ) attributed the discrepancies between 
the data and a pure 3:2 ratio to the action of an additional 
ad-hoc force which moves the frequencies away from their 
'natural' ratio, along a correlation such as the observed one. 

While the correlation anal yzed in BMH 5 treat ed the 
sample of sources as a whole, IZhang et al. (2006a, here- 
after Z06) analyzed the sources separately and showed that 
there are differ ences between sources, and deviations from a 
linear relation. Hb ramowicz et al. | i|2005al lbl hereafter A05a 
and A05b) presented an analysis of the same correlations. 
Fitting a number of sources separately with a linear func- 
tion Uhigh = o,vi ow + b, they obtain a set of a, b pairs (one 
pair per source). The pairs a,b themselves follow a linear re- 
lation b = Aa + B with intercept B close to 1.5. A05a and 
A05b interpret this fact as a property in accordanc e to the 
resonance model of lAbramowicz fc Kluzniakl (|200ll ). 

Here we show that the analysis of separate sources 
presented by A05a and A05b does not add new informa- 
tion to that of the global correlation. We also re-examine 
the frequency correlation with the addition of new data 
from the literature. In particular, the newly discovered 
high- frequency oscillations in Cir X-l l|Boutloukos et all 
120061 ). which extend the correlation to lower frequencies, 
are important to distinguish between the different classes 
of models. We show that the frequencies of the kHz QPOs 
in Cir X-l are inconsistent w ith the resonance model of 
lAbramowicz fc Kluzniakl (120011 ). 



2 ADDITIONAL DATA 

Since the analysis by BMH05, new pairs of kHz QPOs fre- 
quencies have been published, some of which are particu- 
larly important for checking and extending the correlation 
between the kHz QPO frequencies. 



2.1 Cir X-l 

Recently, twin kHz QP Os have been discovered in the pe- 
culiar LMXB Cir X-l (jBoutloukos et all I2006T ). At vari- 
ance with what is observed in Z and atoll sources (see e.g., 
BMH05), the kHz QPO frequencies observed in this system 
are substantially lower than in the Z and atoll sources (in 
the range 56-225 Hz and 230-500 Hz for the lower and upper 
kHz QPO, respectively) and their difference increases with 
frequency. Adding the 1 1 pairs of kHz QPO frequencies from 
iBoutloukos et aL (|2006T ) to the frequency-frequency correla- 
tion of BMH05, we find that, compared to the other kHz 
QPO sources, they follow a very different correlation (see 
Figure [T|, with linear best-fitting parameters a — 2.64±0.26 
and b = 73±28 (excluding the highest frequency point, 
which clearly would tilt a linear relation). We note that 
this relation inte r sects the 3/2 line at v\ — — 64± 29 Hz. 
IBoutloukos et all (|2006T ) mention that for Cir X-l a con- 



stant vi OVJ :v upp — 1:3 ratio is possible. A linear relation with 
a = 3 and b = (dotted line in Figure [TJ fits the data statis- 
tically worse than the b 7^ one (reduced x 2 of 2.4 and 1.3 
respectively); a fit with b — and a free, with a reduced x 2 



of 1.3, yields a = 3.31 ± 0.10, i.e. inconsistent with a value 
of 3. 



2.2 XTE J1807-294 

Another important source to include in the sample is the ac- 
creting millisecond pulsar XTE J 1807 -294, for which pair s 
of kHz QPOs have been discovered bv lLinares et alj (2005). 
It should be noted that some of the peaks found in this 
source are rather broad. However, in seven cases out of 
eight, iLinares et al. (2005) could firmly identify two peaks 
as lower- and upper-kHz QPOs based on correlations with 
low-frequency features. We consider only these seven pairs 
as bona fide kHz QPOs. These points are included in the 
samples analyzed by A05b. Similar values were reported by 
IZhang et all (|2006bt) using the same data. 



2.3 SAX J1808. 4-3658 

For thi s system, the first dis covered accreting millisecond 
pulsar, IWiinands et alj (|2003h found a pair of kHz QPOs 
at ~ 500 Hz and ~ 6 95 Hz, respectively. In addition, 
Ivan Straaten et "all (|2004T ) reported pairs of broad features 
at relatively low frequencies that they interpreted as the 
lower kHz QPO in this system; however, they also mention 
the possibility that these QPOs are a different component 
altogether. We do not include these points in our sample. 
The point corresponding to SAX J1808. 4-3658 is shown in 
Figure [T] as a six-pointed star. 



3 PROPERTIES OF THE GLOBAL 

CORRELATION VS. SINGLE-SOURCE 
CORRELATIONS 

In the following we discuss the significance of correlations 
between linear parameters derived from separate sources. 
We first show the results of simulations based on the existing 
data. This more intuitive approach is then followed by its 
analytical interpretation. 

3.1 The BMH05 sample and single sources 

As a starting point, we take the sample originally published 
by BMH05 and partially used also by Z06, A05a and A05b, 
plus the recent frequencies fr om the millisecond accreting X- 
ray pulsar XTE J1807-294 (ILinares et al.ll2005l . see above). 
Different classes of sources (Atoll and Z) can follow different 
correlations, as shown by BMH05, A05a, and A05b, and the 
same applies to individual sources (see also Z06). However, 
the fact that the overall BMH05 sample of points shows a 
good degree of correlation (we will assume for the remainder 
of the paper that the correlation is linear, but the precise 
functional shape is not relevant for this discussion), indicates 
that the differences between source classes and individual 
sources cannot be very large. 

In Figure [2 we plot the BMH05 sample together with 
the best-fit models found in A05a and A05b for selected 
sources. Notice that discrepant a and b values are quoted 
in A05a and A05b for the same sources; for the discrepant 
cases, we use the A05a values. We remark that Figure[2]con- 
tains also points from sources not fitted with the model lines 
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Figure 1. Plot of the kHz QPO frequency- frequency relation including the points f rom Cir X-l |B outlou kos et alj|2006h . The XTE 
J1807-294 data are from lLinares et al.l J2005l> and the SAX J1808.4-3658 point is from lWiinands et al.l II2003T) . The~asterisks are pairs of 
confirmed simultaneous high-frequency QPOs in black-hole systems. The thin solid line represents a fixed 3:2 ratio, and the dotted line 
a 3:1 ratio. The dashed line is a fit to the Cir X-l points, excluding the highest frequency one. The thick solid li ne shows the relation 
between the periastron-precession and the Keplerian frequencies for M = 2Mq; in the relativistic precession model l|Stella fc Vietrilll999h 
those two frequencies have been identified with the lower and upper kHz QPO, respectively. The dark area represents i>[ on , > Vhigh> 
obviously not allowed by definition. 



plotted in that figure. Unless there is an unlikely correlation 
between dynamic range in frequencies and deviation from 
the average correlation, the moderate spread observed in 
the data excludes large deviations in correlation parameters 
for single sources. In other words, in Fig. [2] the spread of 
the single-source correlations, when extrapolated over the 
full range of the sample, is larger than the spread of the 
points. Slopes very different from that of the sample cor- 
respond to sources for which a reduced dynamic range of 
frequencies is observed. This leads to the conclusion that it 
is likely that these different correlations are the result, at 
least in part, of the limited range of kHz QPO frequencies 
covered by these sources. Notice that, at variance with what 
reported by A05a, the best fit models to different sources do 
not intersect at the same point. 

3.2 The sample as a whole 

As shown by BMH05, the lower and upper kHz QPO fre- 
quencies of a large sample of sources are roughly linearly 
correlated as v upp = 0.92i/; otl) + 360. The precise values are 



irrelevant here, as well as the precise functional shape of 
the correlation (see also Z06) . We consider here the BMH05 
sample without taking into account the association of points 
with separate sources. We simulate artificial sources by tak- 
ing 1000 random sub-samples of n real points where n is also 
a random number for each artificial source, chosen between 
20 and 80. For each sub-sample, we fit the n values with 
a linear function and thus derive a set of 1000 pairs of a 
and b. A plot of these pairs shows a marked anti-correlation 
b = Aa + B with slope A = (-1.69 ± 0.01) x 10~ 3 and in- 
tercept B — 1.51T0.01, consistent with the values quoted 
by A05b from their source-by-source analysis (see Figure 
[3] the range in a and b values in our simulation is limited 
because the procedure does not likely produce random sub- 
samples spanning a limited range in frequency, as it appears 
to happen in reality.) From this simulation we find that any 
randomly selected subset of points from this dataset yields 
a, b pairs that are distributed along such a straight line, and 
have B ~ 1.5. Simulating random straight lines in the plane 
would of course not lead to any a vs. b correlation, but in 
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Figure 2. Frequency-frequency correlation for the sources in 
BMH05, with the addition of the centroid frequencies of XTE 
J1807-294 from Linares et al. (2005) [crosses, see text]. The lines 
are the best fit lines from A05a and A05b. 

this case we are extracting lines that cross the correlation 
segment visible in Figure [2] 

The origin of the correlation shown in Figure [3] can be 
seen more clearly through a simulation that does not make 
use of the actual observed frequencies. We simulated an arti- 
ficial set of pairs of frequencies following a linear correlation 
y — 0.92x + 360 over the range fi = 150-900 Hz and added 
a Gaussian spread with <r=50 Hz. Then, we repeated 100 
times the procedure outlined above (extract 1000 random 
subsamples, fit them with a linear function, fit linearly the 
a values vs. the b values and obtaining the intercept B), 
obtaining a set of 100 values for B. The resulting distribu- 
tion of the 100 B values is strongly peaked around 1.5; this 
shows that, while this value could originate from the differ- 
ences in the frequency-frequency distributions for different 
sources, as suggested by A05a,b, it could also simply be the 
consequence of the values a— 0.92 and &=360 of the global 
correlation and of the ranges spanned by the observed fre- 
quencies. We repeated the simulation with the same a and 
b values, but with v\ in the ranges 600-1000 Hz and 10-500 
Hz, and obtained average B values of 1.36 and 2.0 respec- 
tively. 

3.3 Significance of the anti-correlation between 
linear parameters 

The results of the simulations shown above can be inter- 
preted analytically. The a,b parameters of a y = ax + b 
correlation are covariant. This means that the confidence 
contours around the best fit value are always elongated di- 
agonally in the direction of anti-correlation. The inclination 
of the confidence ellipse is determined both by the best fit 
parameters and the distribution of the points along the cor- 
relation. In particular, the intercept A of the major axis of 
the ellipse can be found simply by fitting a linear relation of 
the form y = Ax to the data pairs Xj , yj . Analytically, this 
can be calculated as 

A y!L±y± (i) 



If the errors are all the same size, Eqn. 1 reduces to 

A =< yx > I < x 2 > (2) 

For the typical range of x — v\ and y = V2 in the data, this 
can be approximated as A =< v-z > / < v\ > Obviously, 
a fit with zero intercept to data linearly distributed with 
non-zero intercept will cross the data around the middle of 
the points, which is roughly centered at < V2 > / < v\ >. 
This means that, given our datasets, the values of A will 
depend also crucially on the distribution of points (see also 
BMH05). 

We selected a clean sample of frequencies from a few 
sources, fitted with a linear function with zero intercept and 
computed the A values from Eqn. 1, which resulted compat- 
ible with the fitting parameters (see Table [TJ where A values 
obtained with Eqn. 1 are shown, as well as < V2 > / < v\ > 
values). Notice that for GX 17+2 and GX 5-1 we only se- 
lected pairs of frequencies where both detections were above 
3cr in Homan et al. (2002) and Jonker et al. (2002a). For 
GX 340+0, this selection from the data in Jonker et al. 
(2000) resulted in only three pairs of frequencies: we there- 
fore excluded this source from the sample. It is clear that the 
A (excluding Cir X-l) values cluster around the value 1.5. 
Fig. [4] shows the best-fitting A and B values and confidence 
contours for these sources. The alignment of all ellipses is 
evident. 

As we have shown above, the values of A cluster around 
the value 1.5 because the observed linearly-correlated fre- 
quencies show a ratio of mean values close to 1.5. This 
was also shown in BMH05. This means that the confidence 
contours will on average have an inclination pointing at 
B = 0,A = 1.5. 



4 DISCUSSION 

Aside from precise fitting of separate sources with differ- 
ent models, it is evident from Figure [2] that the centroid 
frequencies of kHz QPOs of both atoll and Z sources are 
roughly distributed, although with some scatter, around a 
line approximated by Vhigh = 0.92ui om + 360 (see BMH05). 
We have shown that this fact alone leads to a strong anti- 
correlation, with intercept very close to the value 1.5 (the 
exact value of the intercept depends on the range of fre- 
quencies spanned by the QPOs), between the parameters a 
and b obtained from linear fits to any sub-sample of points, 
or to points generated randomly along the global correla- 
tion. This result alone implies that the analysis of separate 
sources cannot be used in the way proposed by A05a,b for 
testing the resonance model. 

The question however remains why the frequencies are 
distributed along that major correlation and with these fre- 
quency boundaries. As shown above, the intercept of 1.5 in 
the a-b anti-correlation could be different should the linear 
parameters of the sample of sources, or the range spanned 
by the kHz QPO frequencies be different. Moreover, notice 
that there are different combinations of a, b and frequency 
range that yield to B ~1.5. The global correlation visible 
in Figure [1] and its frequency range, constitute therefore 
one of the major observational results that must be directly 
addressed by theoretical models. 

Figure [1] shows that the frequency points from Cir X- 
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Name a b A < v% > / < v\ > Ref. 



GX 5-1 


0.856 ± 0.022 


374 ± 8 


1.85 


1.95 


1 


Sco X-l 


0.791 ± 0.005 


430 ± 3 


1.41 


1.42 


2 


GX 17+2 


0.870 ± 0.041 


361 ± 26 


1.45 


1.46 


3 


4U 0614+09 


1.078 ± 0.020 


269 ± 12 


1.49 


1.52 


4,5 


4U 1728-34 


0.899 ± 0.036 


418 ± 25 


1.47 


1.46 


6 


4U 1636-53 


0.672 ± 0.025 


541 ± 22 


1.32 


1.34 


7.8 


4U 1608-52 


0.810 ± 0.027 


419 ± 16 


1.45 


1.42 


6 


4U 1820-30 


0.892 ± 0.045 


355 ± 33 


1.39 


1.39 


9 


4U 1915-05 


1.157 ± 0.045 


256 ± 20 


1.59 


1.62 


10 


XTE J1807-294 


1.145 ± 0.124 


171 ± 30 


1.90 


1.95 


11 


Cir X-l 


2.339 ± 0.474 


104 ± 58 


3.11 


3.18 


12 



Table 1. Sample fo LMXBs with best fit linear parameters, A values obtained with Eqn. 1, < kj > / < fi >, and references. 1. Jonker 
et al. (2002a); 2. Mendez & van der Klis (2000); 3. Homan et al. (2002); 4. van Straaten et al. (2000); 5. van Straaten et al. (2002); 6. 
Mendez et al. (2001); 7. Di Salvo et al. (2003); 8. Jonker et al. (2002b); 9. Zhang et al. (1998); 10. Boirin et al. (2000); 11. Linares et al. 
(2005); 12. Boutloukos et al. (2006). 
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b values 



Figure 3. Linear parameters a and b obtained from 1000 ran- 
dom samples from the points in Figure [2] The points lie on a cor- 
relation which is itself linear with intercept S=1.51±0.01. The 
squares correspond to the a,b pair obtained for Cir X-l including 
the highest-frequency point (full error bars) and without (dashed 
error bars). 



1 pose a more serious pro blem to the resonance model 
of I Abramowicz et al. | l|2003l ): They are distributed along a 
steeper line than the one followed by the frequencies of 
the sources in BMH05, A05a, and A05b, constituting a 
low-frequency extension of the correlation defined by those 
sources. Notice that the linear parameters a and b for Cir 
X-l are very different from those of the other sources, in a 
way that is not at all consistent with the correlation shown 
in Figure 2 of A05b (see Figure [3} . 

It is clear from Figure [T] that the frequencies of Z and 
atoll sources do not follow a re lation Vhi a h = l -5vi ow , as al- 
ready pointed out by BMH05. iRebuscd l|2004l) presented a 
model where the addition of an ad hoc force on an oscillator 
responsible for an otherwise fixed pair of frequencies could 
lead to the observed co rrelation. The kHz QPOs in Cir X-l 
( Boutl oukos et al ]|2006l ) follow a correlation which intersects 




Figure 4. Confidence contours (68.3%, 90%, 99%) corresponding 
to the linear fits to the sources in Table [T] with the exception of 
XTE J1807-294 and Cir X-l. The dashed line is the linear fit 
to the data, yielding a slope of 1.55+0.07 and an intercept of 
-0.0018+0.0002. 



the 3:2 fixed-ratio line at v ~ —64 Hz, i.e. the two lines do 
not intersect in the physical range of interest for the frequen- 
cies. This, together with the fact that the kHz QPOs in Cir 
X-l show other properties and correlations very similar to 
the kHz QP Os in other neut ron-star systems, cast doubt on 
the model of lRebuscd l|2004h . The possibility of Cir X-l be- 
ing a black hole is very unlikely. Eight X-ray bursts were ob- 
served with EXOSAT (jTennant. Fabian fc Shafe3ll986al lbh . 
the high-frequency QPOs are consistent with being neutron- 
star kHz QPOs, and at high luminosities the tracks in the 
X-ray color-color diagrams and the low-fr equency variabilit; 
are v ery similar to those for the Z sources i|Shirev et al.lll99"r 
1999). Also, it is suggestive that the Cir X-l frequencies, 
while following a much steeper linear relation, join smoothly 
with the low-frequency end of the Z sources. 

One of t h e m ost interesting results from 
iBoutloukos et aLI ([2006 ) is that the Cir X-l points are 
the first kHz QPOs observed with a frequency difference 
which increases with frequency, in accordan ce with the pre- 
diction of the relativistic precession model (S tella fc Vietril 
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19991) and th e Alfven wave model l|Zhandl20o3 ; IZhang etafl 
2005, l2006al ). In Fi gure[TJ we also plot the line for the rela- 
tivistic precession model corresponding to a mass of 2M@ 
(the line for the Alfven wave model lies very close to that of 
the relativistic precession model). The agreement with Cir 
X-l and the Z sources at low frequencies is remarkable. At 
high frequencies, the points lie above the model line; this of 
course could reflect different masses for the neutron stars, 
but for som e sources additional assumptions are required 
anyway (see IStella fe Vietrilll999h . Moreover, above i/i = 150 
Hz the data follow a rather straight correlation, while the 
model line is significantly curved at all frequencies. 

The points corresponding to the two accreting mil- 
lisecond pulsars, XTE J 1807-294 and SAX J1808.4-3658 
(|Linares et alj 120051 ; IWiinands et al.1 120031 ) lie below the 
global correlation. It is interesting to note that from their 
analysis of correlations between all timing frequencies of a 
sample of atoll sources and accreting millisecond pulsars 
(not including XTE J1807-294), ? conclude that for two of 
the millisecond pulsars (including SAX J1808. 4-3658) both 
the lower and upper kHz QPOs appear shifted down by a 
factor ~ 1.5 (?, have to multiply V\ ow and v upp of these two 
millisecond pulsars by ~ 1.5 for them to match the corre- 
lation of the other sources in their sample). Applying a 1.5 
factor to both kHz QPO frequencies corresponding to the 
pulsars in our sample, the two sources are indeed moved 
to the global correlation. Within the relativistic precession 
model, such a shift could imply a neutron-star mass higher 
by a factor of 1.5, which is unlikely as it would lead to ac- 
creting millisecond pulsars containing neutron stars of 3Mq , 
assuming 2M© for the other sources. We note that a simi- 
lar shift does not move Cir X-l into the linear correlation 
defined by the other sources. 

On Figure [T] we also plot the four points correspond- 
ing to the pairs of confirmed simultaneous high-frequency 
QPOs from black-hole systems (GRS 1915+105: Strohmayer 
2001a; GRO J1655-40: Strohmayer 2001b; XTE J1550-564: 
Miller et al. 2001; H 1743-322: Homan et al. 2005). Three 
of these points lie very close to the 3:2 line, while GRS 
1915+105 is consistent with a 5:3 ratio (it appears close 
to the 3:2 line due to the low frequencies involved). Clearly, 
for black hole binaries, a model in which the QPO frequency 
ratio is the ratio of two integer numbers is consistent with 
the available data. 



5 CONCLUSIONS 

We showed that the correlation between parameters ob- 
tained from linear fitting to separate sources applied to atoll 
and Z sources by A05a and A05b is a direct result of the ob- 
served general correlation and cannot be used for testing 
their models. Moreover, the newly discovered kHz QPOs in 
the peculiar neutron-star binary Cir X-l follow a different 
correlation than atoll and Z sources. This observed correla- 
tion appears to be in contradiction with the predictions of 
all modifications of resonance models, but is in good general 
agreement with the relativistic precession model (although 
it is known that the model in its simplest version does not 
explain all observables, such as the relation between kHz 
QPOs and spin frequency, and needs correction at high fre- 
quencies) and the Alfven wave model. Our results indicate 



that the observed high frequency QPOs in black hole and 
neutron star low-mass X-ray binaries most likely have a dif- 
ferent origin. 
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